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UNIFORM REFRACTION IN NEGATIVE REFRACTIVE INDEX 

MATERIALS 

CRISTIAN E. GUTIERREZ AND ERIC STACHURA 


Abstract. We study the problem of constructing an optical surface separating 
two homogeneous, isotropic media, one of which has a negative refractive index. 
In doing so, we develop a vector form of Snell's law, which is used to study 
surfaces possessing a certain uniform refraction property, both in the near and 
far field cases. In the near field problem, unlike the case when both materials 
have positive refractive index, we show that the resulting surfaces can be neither 
convex nor concave. 
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1. Introduction 

Given a point O inside medium I and a point P inside medium II, the question 
we consider is whether we can find an interface surface T separating medium I 
and medium II that refracts all rays emanating from the point O into the point 
P. Suppose that medium I has index of refraction n.\ and medium II has index 
of refraction ft 2 - In the case that both n.\, n .2 > 0, this question has been studied 
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extensively (see in particular 0] and the references therein). This is the so called 
near field problem. The far field problem is as follows: given a direction m fixed 
in the unit sphere in IR 3 , if rays of light emanate from the origin inside medium I, 
what is the surface T, interface between media I and II, that refracts all these rays 
into rays parallel to ml In the case that both n.\, n 2 > 0 this problem has also been 
studied extensively Q. 

A natural question in this direction is what happens in the case when one of the 
media has a negative index of refraction? These media, as of now still unknown to 
exist naturally, constitute the so called Negative Refractive Index Materials (NIM) 
or left-handed materials. The theory behind these materials was developed by 
V. G. Veselago in the late 1960's Hl3l , yet was put on hold for more than 30 years 
until the existence of such materials was shown by a group at the University of 
California at San Diego Ifllll . 

This note is devoted to studying surfaces having the uniform refraction property 
(as discussed above) in the case when k := < 0, both in the far field and near 

field cases. We begin discussing in Section 2.1 the notion of negative refraction, 
and in Section |Z2| we formulate Snell's law in vector form for NIM. Next, in Section 
[3] we find the surfaces having the uniform refraction property in the setting of 
NIMs. In contrast with standard materials, i.e. those having positive refractive 
indices, the surfaces for NIMs in the near field case can neither be convex nor 
concave; this is analyzed in detail in Section |3~4} 

Finally, in Section [4] using the Snell law described in Section |2.2[ we study 
the Fresnel formulas of geometric optics, which gives the amount of radiation 
transmitted and reflected through a surface separating two homogeneous and 
isotropic media, one of them being a NIM medium. We expect that the study in 
this paper will be useful in the design of surfaces refracting energy in a prescribed 
manner for NIMs as it is done for standard materials in 0], 0, and (6]|. 


2. Geometric Optics in Negative Refractive Index Materials 

2.1. Introduction to Negative Refractive Index Materials. The notion of negative 
refraction goes back to the work of V. G. Veselago in Ifl3l . In the Maxwell system of 
electromagnetism, the material parameters e, p are characteristic quantities which 
determine the propagation of light in matter. If E = A cos(r • k + cot), then in the 
case of isotropic material the dispersion relation (7[ Formula (7.9)] is given by 



with n 2 the square of the index of refraction of the material, which takes the form 


(2.1) n 2 = ejj. 

Veselago showed that in the case when both e, p < 0, we are forced by the Maxwell 
system to take the negative square root for the refractive index, i.e. 
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(2.2) n - - yjey. 

In particular he observed that a slab of material with 

(2.3) e = -1 and p = -1 

would have refractive index n = -1 and behave like a lens. A remarkable prop¬ 
erty of this material was shown in 0: that the focusing is perfect provided the 
condition (|2.3|) is exactly met. The theory behind NIMs opened the door to study 
the so called metamaterials, see e.g. IfTOl and ||2|]. 

We end this section by making some remarks on some important differences in 
geometric optics in the case of NIMs. One must generalize the classical Fermat 
principle of optics, see e.g. jit Sec. 3.3.2], to handle materials with negative 
refractive index. This was done not too long ago by Veselago in ffl4l . In particular, 
in the setting of NIMs, the optical length of a light ray propagating from one 
NIM to another is negative because the wave vector is opposite to the direction 
of travel of the ray. However, in this case, it is not possible a priori to determine 
that the path of light is a maximum or minimum of the optical length, as this 
depends heavily on the geometry of the problem. Finally, the Fresnel formulas of 
geometric optics must also be appropriately adjusted for negative refraction; we 
analyze this in detail at the end. 

2.2. Snell's law in vector form for k < 0. In order to explain the Snell law of 
refraction for media with k < 0, we first review this for media with positive 
refractive indices. 

Suppose T is a surface in 1R 3 that separates two media I and II that are homoge¬ 
neous and isotropic, with refractive indices Hi and n 2 respectively. If a ray of lightQ 
having direction x G S 2 , the unit sphere in IR 3 , and traveling through medium I 
strikes I at the point P, then this ray is refracted in the direction m G S 2 through 
medium II according to the Snell law in vector form: 

(2.4) ri\{x x v) = n 2 (m x v), 

where v is the unit normal to the surface to f at P pointing towards medium II; 
see |[8l Subsection 4.1], It is assumed here that x ■ v > 0. 

This has several consequences: 

(a) the vectors x, m, v are all on the same plane (called the plane of incidence); 

(b) the well known Snell's law in scalar form holds: 

Hi sin 0i = n 2 sin 0 2 , 

where 0i is the angle between x and v (the angle of incidence), and 0 2 is 
the angle between m and v (the angle of refraction). 

'Since the refraction angle depends on the frequency of the radiation, we assume that light rays 
are monochromatic. 
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Equation ( 2A) is equivalent to (;qx - n 2 m) x v = 0, which means that the vector 
H\X - n 2 m is parallel to the normal vector v. If we set k = h 2 /hi, then 

(2.5) x — k m = Av, 


for some A G 1R. Notice that (2.51 univocally determines A. Taking dot products 
with x and m in ( |2.5| > we get A = cos 6\ - k cos G 2 , cos G\ = x ■ v > 0, and cos G 2 = 
= -\/l - k _ 2 [1 — (x • v) 2 ]- In fact, there holds 


m ■ v 


( 2 . 6 ) 


A = x -v 


K 


Vf 


K 


-2 


(1 - (x • v) 2 )- 


It turns out that refraction behaves differently for 0 < k < 1 and for k > 1. 
Indeed, 

(1) if 0 < k < 1 then for refraction to occur we need x ■ v > Vl - k 2 , and in (2.6 > 
we have A > 0, and the refracted vector m is so that x ■ m > k; 

(2) if k > 1 then refraction always occurs and we have x ■ m > 1/k, and in (2.6) 
we have A < 0; 

see |0 Subsection 2.1]. 

We now consider the case when either medium I or medium II has negative 


refractive index, so that k 


ni 

Hi 


< 0. Let us take the formulation of the Snell law 


as in (2.51. This is again equivalent to 
(2.7) x-Kin 


Av, 


with x, m unit vectors, v the unit normal to the interface pointing towards medium 
II. We will show that 


( 2 . 8 ) 


A = x-v + 


In fact, taking dot products in (2.7), first with x and then with m, yields 
l-KX-m = Ax-v, and x ■ m — k - Am - v. 


1 1 

This implies that- Ax ■ v - x • m 

K K 


k + Am ■ v. We seek to get rid of the term 


m ■ v. To do this, from (2.71 we have that m 


x - Av 


K 


so m ■ v 


substitution, we therefore obtain the following quadratic equation in A: 

A 2 - 2Ax • v + (1 - k 2 ) = 0. 


1 1, „ 

—x ■ v-A. By 

K K 


Solving this equation yields 

A = X -V ± yj{x ■ v) 1 - (1 - K 1 ). 

If —oo < k < —1, then (x • v) 2 - (1 - k 2 ) > 0 and the square root is real. On the other 
hand, if — 1 < k < 0, then the square root is real only if x • v > Vl - k 2 which means 
that the incident direction x must satisfy this condition. It remains to check which 
sign (+) to take for A. Recalling that 


x-v - Km ■ V - A - X ■ V ± yj(x ■ v) 2 - (1 - K 2 ), 
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we see, since k < 0 and m ■ v > 0, that we must take the positive square root. 
Hence we conclude (2.8). Notice that in contrast with the case when k > 0, the 
value of A given in (2.8 1 is always positive when k < 0, and it can be also written 
in the following from 


A 


x • v 


+ |k| a/1 - K ~ 2 (1 - (x ■ v) 2 ). 


which is similar to (2.6). In the "reflection" case k = -1, since x-v > 0, this formula 
yields 

A = 2 x • v, 

and hence the "reflected" vector m is given by 

m = 2 (x • v) v - x. 

That is, after striking the interface T the wave with direction m travels in the 
material with refractive index n 2 . 

Note that the vector formulation ( |2.7| ) is compatible with the Snell law for 
negative refractive index ( IH5ID . Indeed, taking the cross product in ( 2 .7) with the 
normal v yields x X v = k m X v, and then taking absolute values yields 

sin Q\ = —k sin d 2 , 

where Q\ is the angle between x and v, and d 2 is the angle between m and v. 



Figure 1. Snell's law in the case of normal refraction and negative 
refraction (dotted arrow). 
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2.3. Physical condition for refraction when k < 0. We analyze here conditions 
under which an electromagnetic wave is transmitted from medium I to II and 
there is no total internal reflection. 


2.3.1. Case when —1 < k < 0. The maximum angle of refraction is 0 2 = n/2 which 
is attained when sin Q\ = —k, that is, the critical angle is 0 C = arcsin(-jc). If 6 \ > 0 C , 

there is no refraction. So 0 2 = arcsin |— sin G\ j, for 0 < Q\ < G c . The dot product 
x ■ m - cos(@i + 6 2 ). Let h{Q\) = Q\ + 0 2 = @1 + arcsin sin Q\ j. We have 


fe'(0i) = 1 - - cos G 1 1 

k 


> 0 


sin GJk 2 

for 0 < @1 < Q c , and so h is increasing on [0, 6 C ] and therefore 

G\ + @2 — G c + 71 / 2 , for 0<@!< G c . 


Then the physical constraint for refraction is 

(2.9) x ■ m = cos(@i + @ 2 ) > cos(0 c + n/2) = - sin 6 C = k. 


2.3.2. CasewhenK < -1. InthiscaseO < Q\ < n/2, and since 62 = arcsin sin d\ 

the maximum angle of refraction is when 6 1 = n/ 2 , that is, the maximum an¬ 
gle of refraction is 6 r = arcsin(-l/x). Now h{ 6 \) is increasing on [0, ti/2], so 
0i + 62 < n /2 + arcsin(-l/K), and 

1 

(2.10) x - in = cos(0i + 0 2 ) > cos {n /2 + arcsin(-l/K)) = -. 

K 


3. Surfaces with Uniform Refraction Property 


3.1. The Far Field Problem. Let m e S 2 be fixed. We ask the following question: 
if rays of light emanate from the origin inside medium I, what is the surface T, 
interface between media I and II, that refracts all these rays into rays parallel to ml 
Here I has refractive index n lr II has refractive index n 2 such that k = n 2 /n l < 0. 

This question can be answered using the vector form of Snell's law. Suppose 
T is parameterized by the polar representation p{x)x where p > 0 and x G S 2 . 
Consider a curve on T given by r(t) = p(x(t))x(t) for x(t) G S 2 . According to ( |2.7| >, 
the tangent vector r'(t ) to T satisfies r'(t) ■ (x{t) - Km) = 0. That is, ([p(x(t))]'x(t) + 
p(x(t))x'(t)j ■ (x(t) - Km) = 0, which yields (p(x(t))(l - Km ■ x(t)))' = 0. Therefore 


(3.11) 


P(x) = 


1 - Km-x 

for x G S 2 and for some b G R. To see what is the surface described by ( |3.11[ ), we 
assume first that —1 < k < 0. Since m ■ x > —1, we have 1-Kra-r > 1 +k> 0, and 
so b > 0. 
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Suppose for simplicity that m = e 3 , the third-coordinate vector. If y = ( 1 /, y 3 ) e 
R 3 is a point on T, then y = p(x)x with x = y/|y|. From ( j3.HI ), |y|—7c y 3 = b, that is, 
|y'| 2 +y 2 = (k y 3 + b ) 2 which yields |y'| 2 +(l - x 2 )y 2 - 2 xby 3 = b 2 . This equation can 
be written in the form 


(3.12) 


I y 


/12 


ys 




1 - K 2 , 


+ 


vr 




1 




which is an ellipsoid of revolution about the y 3 -axis with upper focus (0,0), and 
lower focus (0,2x&/(l - k 2 )). Since |y|= rcy 3 + b and the physical constraint for 

y 

refraction (|2.9b, — • e 3 > k is equivalent to i/ 3 > ---. That is, for refraction 

1 — 1 |y| 1 - k 2 

to occur y must be in the upper half of the ellipsoid ( ]3. 12[ > ; we denote this semi¬ 
ellipsoid by E(e 3/ b). To verify that E(e 3 , b) has the uniform refracting property, that 
is, it refracts any ray emanating from the origin in the direction e 3 , we check that 


\2.7) holds at each point. Indeed, if y e E(e 3 , b), then I p- - xe 3 


lyl 


y_ 

lyl 


y 

1 - KC 3 • —r > 

lyl 


xe 3 \ 


y 

e 3 > 0, and so — 

lyl 


xe 3 is an outward normal to E(e 3 , b) 


1 + k > 0, and (- 

\lyl 

at y. 

Rotating the coordinates, it is easy to see that the surface given by p.llj ) with 
-1 < k < 0 is an ellipsoid of revolution about the axis of direction m with upper 

focus (0,0) and lower focus -- -m. Moreover, the semi-ellipsoid E(m,b) given 

1 K 

by 


(3.13) 


E(m,b ) = < p(x)x : p(x) = 


xm- x 


x e S 2 , x-m>K\ 


has the uniform refracting property: any ray emanating from the origin O is 
refracted in the direction m, where — 1 < k < 0. See Figures [2] and [3] comparing 
refraction when k < 0 versus k > 0, where P = (0,2icfo/(l - k 2 )) and m = e 3 . 

Let us now assume that K = -1. From ( 3.11| ) we have |y|+y 3 = b, and so 
|yf+y 2 = b 2 - 2by 3 + y 2 , and therefore 


y3 = ^(t 2 -iyf)- 


This is a paraboloid with axis e 3 and focus at O. Since k = —1, from the physical 
constraint for refraction (2.9) we get any ray with direction X £ S 2 emanating from 
the origin is refracted by the paraboloid. 

Now we turn to the case k < — 1. Due to the physical constraint of refraction 
(2.10|), we must have b > 0 in (3.11); we determine what kind of surface this is. 
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Figure 2. Ellipsoid re¬ 
fracting when -1 < 

k < 0 


Figure 3. Ellipsoid re¬ 
fracting when 0 < k < 
1 


Define for b > 0 

(3.14) H(m,b) = i p(x)x : p(x ) 


1 - Km-x 


, x £ S 2 , x ■ m > 1/k> . 


We claim that H(m, b) is the sheet with opening in direction m of a hyperboloid 
of revolution of two sheets about the axis of direction m. To prove the claim, 
set for simplicity m = e 3 . If y = (y', 1 / 3 ) e H(e^,b), then y = p(x)x with x = 
y/|y|. From ( |3.14| >, -k y 3 + \y\= b, and therefore \y'\ 2 +y\ = (k y 3 + b) 2 which yields 

y}) \ / 

b 2 . Thus, any point y on H(e 3 , b) satisfies 


|y'| 2 -( 7 c 2 -l) 
the equation 

(3.15) 


V 3 + 


k a 


1 


, 2/3 + 


xb 


K 


I y 


f\2 


. K 2 - 1, 


Vk 2 - 1/ 
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which represents a hyperboloid of revolution of two sheets about the y 3 axis with 
foci (0,0) and (0, -2 xb/(x 2 - 1)). The sheets of this hyperboloid of revolution are 
given by 

xb 

ys 


K 2 ~ 1 7C 2 — 1 


1 + 


I y 


/|2 


A (&/ Vk 2 - 1 ) 


We decide which one satisfies the refracting property The sheet with the minus 

b 

sign in front of the square root satisfies xy 2 + b > -- > 0, and hence has 


K 


1 


and therefore this is the sheet to consider. For 


polar equation p(x ) = -- 

-L K. £3 * % 

a general m, by a rotation, we obtain that H(m, b) is the sheet with opening in 

direction opposite to m of a hyperboloid of revolution of two sheets about the axis 

_ 

of direction m with foci (0,0) and —-— -m. 


K z 


1 


Notice that the focus (0,0) is inside the region enclosed by H(m, b) and the focus 

—2 xb y 

m is outside that region. The vector — - Ktn is an outer normal to H(m, b) 


X 


1 


at y, since by ( |3.14| ) 


lyl 


irH 


-2xb 
, k 2 - 1 


-m 


y U 


-2b 2x 2 b 

+ 


K 2 — 1 


-2b 


X 


1 


+ 


K 2 ~ 1 

2 x 2 b 

K 2 - 1 


+ xm - y - |y| 
- b = b > 0. 


y \ 1 y \ y 

Clearly, — - xm ■ m > - x > 0 and — - xm • — > 0. Therefore, H(m, b) 

\lyl / k V lyl /lyl 

satisfies the uniform refraction property. We summarize the uniform refraction 
property below. 

Theorem 3.1. Let n\ r n 2 be two indices of refraction for media I and II, respectively, and 
set x = n 2 /ri]. Assume the origin is in medium I, and E(m, b), H{m, b) defined by ( |3.11[ ) 
and ( |3.14[ ) respectively. Then: 

(1) If—1 < x < 0 and E(m, b) is the interface between medium I and medium II, then 
E(m,b ) refracts all rays from the origin O into rays in medium II with direction 

m. 

(2) If x < -1 and H(m, b) is the interface between medium I and medium II, then 
H(m,b) refracts all rays from the origin O into rays in medium II with direction 

m. 


3.2. The Near Field Problem. Given a point O inside medium I and a point P 
inside medium II, we construct a surface F separating medium I and medium II 
such that T refracts all rays emanating from O into the point P. Suppose O is the 
origin, and let X(t) be a curve on T. Recall Snell's law says that 

x - xm = Av 
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where v is the normal to the surface T, and we are considering the case when 
k < 0. Then via Snell's law, we see that the tangent vector X'(t) must satisfy 


That is to say. 

Hence T is the surface 


X\t) 


(m p - m \ 

\|X(t)| |P-x(t)|/ 

ix(or+Kip-x(t)r=o. 


= 0 . 


(3.16) 


|X|+k|X-P|= b. 


Compare this surface to the Cartesian ovals which occur in the near field case 
when k > 0. In this case, since the function F(X) = |X|+k|X - P| is convex and the 
ovals are the level sets of P, then the ovals are convex. However, in the case when 
k < 0, the function F(X) is no longer convex and its level sets are convex sets only 
for a range of values in the parameter b. That is, the surface ( 3.16 [ > is in general 
not a convex surface, see Theorem |3.4| and Figure [5 

Case — 1 < K < 0. We will show that the surface ( 3.1 6[ ) is nonempty if and only 
if b is bounded from below. In fact, this follows by the triangle inequality; since 
1 + k > 0, we see that |X|+k |X - P|> k |P| for all X. Therefore if the surface ( 3.16 ) 
is nonempty then 

(3.17) b > k\P\. 


Viceversa, if b satisfies ( |3.17| ), then the surface ( |3.16[ ) is nonempty. Let X 0 be of the 
form X 0 = AP. We will find A > 0 so that X 0 is on the surface defined by ( |3.16| ). 
This happens if A \P\+k |A - 1| |P|= b. Suppose b > |P|. Then we must have A > 1. 
So b = A \P\+k (A - 1) |P|, and solving for A gives 

b + k\P\ 


A = 


Now suppose k\P\< b < |P|, and let 


A = 


(1 + K) |Pf 


b - k\P\ 


(1 - K) |P| 

We then have 0 < A < 1 and X 0 = A P is on the surface ( |3.16| >. [] 

Case k = -1. In this case to have a nonempty surface we show that b must be 
bounded above and below. If the surface is nonempty then by application of the 
triangle inequality we get that 

(3.18) - |P|< b < |P|. 

Viceversa, if ( |3.18| ) holds, then the surface is nonempty. In fact, the point X 0 = AP 
. i 1 b + P\ 

with A = _ belongs to the surface. 


2\P\ 


+ Notice that the set {X : |X|+k|X - P|< k|P|+ 6}, 5 > 0 is contained in the the ball with center zero 
and radius 6/(1 + k). 
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Case k < — 1. We will show that the surface ( 3.16[ > is nonempty if and only if 
b is bounded above. First notice that by the triangle inequality and the fact that 
K < — 1 we have that |X|+k |X - P|< |P| for all X, and therefore if (3.16 > is nonempty 
then 


(3.19) 


b < |P|. 


On the other hand, if (|3.19|) holds then will find Xo = AP G F for some A > 0. We 

\) + k\P\ 

want A|P|+k|P||A - 1|= b. If we let A = ——, then A > 1 and the point AP is on 


the curve. @ 


|P|+k|P| 


3.3. Polar equation of (|3.16|). Next we will find the polar equation of the refracting 
surface ( |3.16[ > when k < 0. We show that only a piece of this surface does the 
refracting job, otherwise total internal reflection occurs. This follows from the 


physical conditions described in Section 2.3.1 
Let X = p(x)x with x G S 2 . Then we write 

x\p(x)x - P|= b - p(x). 

Squaring both sides yields the following quadratic equation 
(3.20) (1 - k 2 ) p(x) 2 + 2 (k 2 x ■ P -bj p(x) + b 2 - k 2 |P| 2 = 0. 

Solving for p we obtain 

(b - k 2 x ■ P ) ± yj{b - k 2 x ■ P) 2 - (1 - K 2 )(b 2 - k 2 |P| 2 ) 


p(x) = 


K z 


Define an auxiliary quantity 

A(f) := (b - K 2 t) 2 - (1 - K 2 )(b 2 - k 2 |P| 2 ) 


so that 


P±(x) 


(b - k 2 x ■ P) ± ^A(x ■ P) 

1 - K 2 


3.3.1. Case -1 < k < 0. We begin with a proposition. 

Proposition 3.2. Assume -1 < k < 0. If x G S 2 , then 
(3.21) A(x ■ P ) > k 2 (x ■ P - b) 2 , 

with equality if and only if \x ■ P|= \P\. Hence, the quantity y/A(x ■ P ) is well defined. 


^Notice that the set £ = (X : |X|+k|X - P|> |P|-6}, 6 > 0, is contained in the ball with center 0 
and radius -5/(1 + k). 





















12 


C. E. GUTIERREZ AND E. STACHURA 


Proof. First notice that \x ■ P |< |P| comes for free since |x|= 1. We have, since 
0 < k 2 < 1 and x-P < \P\, 

A(x ■ P ) - k 2 (x ■ P -b) 2 = -2 bx 2 x ■ P + k 4 (x ■ P) 2 + x 2 b 2 

- k 4 |P | 2 + k 2 |P | 2 - k 2 (x ■ P ) 2 + 2k 2 bx ■ P - K 2 b 2 
(3 22) = k\x ■ P) 2 - k 4 |P | 2 + k 2 |P | 2 - k 2 (x ■ P) 2 

= -K 4 (|P | 2 - (X ■ P) 2 ) + K 2 (|P | 2 - (X ■ P) 2 ) 

= K 2 (l - K 2 ) (|P | 2 - (X ■ P) 2 ) 

> 0 

as desired. □ 


Now, notice that in this case, in order to make sense of the physical problem, we 
must have P lying outside of the oval; that is, outside of the region \X\+k\X-P\< b. 
In addition to the requirement that the surface be nonempty, we see that the 
condition on b in order that the ovals are meaningful is that 

k\P\< b < \P\ 


Recalling that 

K\p(x)x - P|= b - p(x) 

we see that, since k < 0, we must have p(x) > b. We now choose which sign (±) to 
take in the definition of p+ so that p > b. In fact, by Proposition [3. 2[ 


P+(x) = 


(b - k 2 x ■ P) + VA(x • P) (b - k 2 x ■ P) + \K\\b -x-P\ 


1 - K 2 


> 


K z 


>b, 


where the last inequality holds since |k|< 1. Therefore the equality p+(x) = b holds 
only if \x • P|= \P\ and b = x ■ P. Also 


(b - k 2 x ■ P) + |k ||b -x-P | 1 j(b - k 2 x ■ P) + \x\(b - x ■ P) for b > x ■ P 

1 - k 2 1 -k 2 \(b - k 2 x ■ P) + |k|(x • P -b) for b < x ■ P 

> b. 


Hence p+(x) > b for all x and b such that b f x-P. From \3.17) the oval is not 
degenerate only for k|P|< b. Reversing the above inequalities one can show that 

p~(x) < b. 

Thus the polar equation of the surface T is given by 

(b - k 2 x ■ P ) + -y/A(x ■ P ) 


(3.23) 


h(x, P, b) = p+(x) = 


K z 


provided that k|P|< b. Furthermore, from the physical constraint for refraction 

_ P - p+(x)x _ 

\2.9) , in this case with in = —-——, so x-m> k, and using the equation ( |3.16| >, 


we must have 


IP - p+(x)x | 

x-P >b. 
















UNIFORM REFRACTION IN NIMS 


13 


n 2 



Figure 4. k - -.7, P = (0,2.5). 


Ending with the case when — 1 < k < 0, given a point P € R 3 and k\P\< b, a 
refracting oval is the set 

(3.24) 0(P, b ) = { h(x, P,b)x:xeS 2 ,x-P>b] 

where 


h(x, P, b) 


(b - k 2 x ■ P ) + y/(b - k 2 x ■ P) 2 - (1 - x 2 )(b 2 - k 2 |P| 2 ) 


1 - K 2 


Figure [4] illustrates an example of such a refracting oval. 


3.3.2. Case k < —1. In this case we must have O not on the oval, that is, O must 
lie outside of the set |X|+k|X - P |> b. Flence in combination with the requirement 
that the surface be nonempty we must have 

k|P|< b < |P|. 


Notice this is the same condition on b as in the case — 1 < k < 0. 
Multiplying ( ]3.20[ ) by -1 and solving for p+ yields 

/ s k 2 x-P -b ± J(b — k 2 x ■ P) 2 - (k 2 - l)(K 2 \P\ 2 -b 2 ) 

- 

To have the square root defined we need 

A(x ■ P) = (b - k 2 x- P) 2 - (k 2 - 1)(tc 2 |P| 2 — & 2 ) > 0. 


From ( 3.19| >, to have a nonempty oval we must have b < |P|. To find the directions 
x for which p is well defined, we need to find the values of t for which 


A (t) = (b - k 2 t) 2 - (k 2 - 1)(k 2 |P| 2 -& 2 ) > 0. 


















14 


C. E. GUTIERREZ AND E. STACHURA 


Let z = b - x 2 t, and then we want t such that 

z 2 > (k 2 - 1)(k 2 |P| 2 -& 2 ). 

If K 2 |P| 2 -b 2 < 0, then this is true for any t. On the other hand, k 2 |P| 2 —b 2 > 0 if and 
only if \b\< -k|P|. So for b with \b\< —k\P\ we have to choose t such that z = b — x 2 t 
satisfied 

|z|> V(k 2 - 1)(k 2 |P| 2 -& 2 ). 

If b - x 2 t > 0, i.e., t < b/x 2 , then b - x 2 t > y/(x 2 - 1)(k 2 |P| 2 -& 2 ), that is, 

b- V(k 2 - 1)(k 2 |P| 2 -& 2 ) 


t < 


k a 


If b - k 2 t < 0, i.e., t > b/x 2 , then x 2 t-b> y/(x 2 - 1)(k 2 |P| 2 -& 2 ), that is, 

b+ V(k 2 - 1)(k 2 |P| 2 -& 2 ) 


(3.25) 


t > 


x A 


From ( 2.10| >, to have refraction in this case, we need to have 

P - xp±(x) 1 


x 


> 


|P - xp±(x)\ X 
Note that this is equivalent to requiring that 

x 2 x ■ P -b > (x 2 - 1 )p±(x) 

That is, we need 


P±(x) < 


x 2 x ■ P -b 
x 2 - 1 


Note that since x 2 > 1, we see that p±(x) < 0 for x 2 x ■ P - b < 0. Thus we must have 


xx 


■ P > b if we are to have p±(x) > 0. From (3.25|> we need that 


x ■ P > 


b+ V(k 2 - 1)(k 2 |P| 2 -& 2 ) 


x 


so that A(x • P) > 0. We are left with choosing either p+(x) or p_(x) now. But by the 
physical restraint for refraction, we see only p-(x) will do the job since 


p-(x) 


x 2 x 


P -b - tJA(x ■ P) x 2 x ■ P — b 


provided that 


x ■ P > 


x 2 - 1 x 2 -1 

b+ V(k 2 - 1)(k 2 |P| 2 -£> 2 ) 


x 


^Notice the requirement |i>|< —k|P| is satisfied since we are assuming that b < |P| so in particular 
I b\< |P|< -k|P|. 
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Hence for k < — 1, refraction only occurs when b < |P| and the refracting piece of 
the oval is given by 


(3.26) 


0(P, b) = \ h(x, P,b)x : x ■ P > 


b+ 


K z 


with 


7 ./,. Tt U\ — „ /,A _ ^ • p - b ) - V( k2 * • P-b) 2 - (x 2 - l)(K 2 \P\ 2 -b 2 ) 
h(x, P, b) 0 —(x') z - . 

k z - 1 

A picture of the refracting piece of this oval can be obtained from Figure [4] by 
reversing the roles of P and O and changing the direction of the rays; see also the 
explanation at the end of the proof of Theorem |3.4[ 


3.3.3. Case x = -1. In this case we have | p(x)x - P|= p(x) - b, so squaring and 
solving for p yields 

b 2 - \P\ 2 


p{x) 


2 {b-x-Py 


To have a nonempty surface, b must satisfy ( |3.18[ ). Since p > 0 and the numerator 

in p is nonpositive, we must have b < x ■ P. The last inequality follows from the 

_ P - p(x)x 

physical condition for refraction \2.9) with k = -1 and m = 


:• Notice that 


p(x) > b, that is. 


|P | 2 


IP - p(x)x\ 

> b, since this is equivalentto(&-x-P) 2 +|P| 2 -(x-P ) 2 > 0. 


2 (b - x ■ P) 

Therefore when k = -1 the surface refracting the origin O into P is given by 


(3.27) 


£{P,b) = 


b 2 - |P | 2 
2 (b-x-P) 


x : x £ S 2 , x ■ P > b\ 


with \b\< |P|. 

The uniform refraction property is summarized below. 


Theorem 3.3. Let ri \, n 2 be two indices of refraction of two media I and II, respectively, 
such that k := n 2 /ni < 0. Assume that O is a point inside medium I, P is a point inside 
medium II, and b £ (k|P|, |P|). Then, 


(1) If—I < k < 0 and I := &(P, b ) is given by [3.24), then T refracts all rays emitted 
from O into P. 

(2) Ifx < -1 and Y := f7(P, b) is given by ( |3.26[ ), then T refracts all rays emitted fro 
O into P. 

(3) Ifx = -1 and Y := S(P, b) is given by ( | 3.27 ), then T refracts all rays emitted fro 
O into P. 


m 


m 
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3.4. Analysis of the convexity of the surfaces. We assume k < 0 and analyze the 
convexity of the surface 

|X|+k|X-P|=&. 

We will show the surface is convex for a range of b's and neither convex nor 
concave for the remaining range of b's. |^] We begin by letting X = {x\, ,x n ) and 
P - (Pi> ’ ’ • ,p n ) and analyzing the function 


F(X) = |X|+k|X-P|. 


We have the gradient 


. . X X - P 
DF(X) = — + k p^Tpj/ 


and 


r _ * 1 XiX J , „ x 1 (*i ~ Pi)( x j ~ Pj) 

XiX < ~ ° ij |X| |X| 3 + K dii |X - P\ K |X - P| 3 


where t) /; is the Kronecker delta. Then 

n 

Q(x,« = 


hi =1 

1 


_ 

|X| |XP 


u 


+ K 


l^p 


K 


bj =1 


|X-P| |X-P| 3 ^ 


u 

- Pi)(Xj - Pj)Z£j 


!i!!_ J_, Y r\2 , , l^l 2 
|X| |X|3 (X ^ 

1 


K 


|X-P| |X-P| 3 


;«X-P)-Z ) 2 


= K 


\X - P I 3 


[Rl 2 IX - P| 2 -((X - P) • a 2 ] + [lil 2 |x| 2 -(x • 4) 2 ], 


By Cauchy-Schwartz's inequality \v ■ w\< \v\ \w\ with equality if and only if v 
is a multiple of zv. If we set £ = X in the quadratic form we get Q(X, X) = 

k [|X| 2 |X - P| 2 -((X - P) • X) 2 ]. Since P / 0, X is not a multiple of X - P 

|X - P | 3 

if and only if X is not a multiple of P and in this case we get Q(X,X) < 0 
since k < 0. On the other hand, if we set E, = X — P, then Q(X, X - P) = 
1 

[|X - P| 2 |X| 2 —(X • X - P) 2 ], and again if X is not a multiple of P we obtain 

Q(X, X - P) > 0. Therefore the quadratic form Q(X, £) is indefinite and therefore 
the function P is neither concave nor convex at each X which is not a multiple of 
P. 

Let now X have the form X = AP, 0 <A<oo, and see for what values of A 
there holds that Q(AP, ^) > 0 for all £. In other words, we study which values of A 
make the Hessian F XjXj (AP) > 0. In particular we wish to find the range ofb that give a 


'To be precise, by convexity we mean that the set that is enclosed by F(X) = |X|+k|X - P|= b is 
a convex set. 
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Figure 5. Ovals k — —.5 


convex surface. 
we have 


The range we obtain will be proven rigorously 


in Theorem 3.4 First off 


Q(AP,£) 


= K iuujw t (1 ~ A)2|4|2|p|2 - (1 - A)2(p ■ « 2 ] + jw ( a2 ' 4 |2|p|2 - a2(p - « 2 1 
= (x +K ^)^[ |tf|P|2 - (P '« 2 l"° 


if and only if 

« A>: =i +K T^ aa 

Let us assume-1 < k < 0. Then^)(A) > 0on(0, 1 /( 1 - k )), <p(A) < 0 on (1/(1-rc), 1), 
cp(A) < 0 on (1,1/(1 + k )), and cp(A) > 0 on (1/(1 + k ), +oo). Therefore 


Q(AP,£) 



for A G (0,1/(1 - k )) U (1/(1 + k ), +oo) 
for A G (1/(1 - k), 1) U (1,1/(1 + k)). 
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Suppose b > k\P\, see ( 3.17) . If A e (0,1/(1 - k)), then since k < 0, we have 
A < 1. If the point X = AP is on the surface ( |3.16| > / then A|P|+k (1 - A)|P|= b and so 

^ ipi 

A = —-rr^r- Therefore when 0 < A < 1/(1 - k) to have F XjXj (AP) > 0, we must 


have 


(l-K)|Pf 


b — k\P\ 1 


(1 - k)\P\ ~ l-K 

That is, 

b < (1 + k)\P\. 

So F XiXj (AP) > 0, 0 < A < 1/(1 - k), with AP on the surface A|P|+k(1 - A)|P|= b if 
and only if b e [k|P|, (1 + k)|P|]. 

Suppose now A G (1/(1 - k), 1). Let us see for which b we have F XiXj (AP) < 0. 

|p| 

we see that in order to have F XiXj (AP) < 0, we need 


Since A = 


(1 - K)|P| 


1 & - k|P| „ 

< T7-!t=7 < 1 


l-K (1 - k)|P| 

so that 

(1 + k)|P|< b < \P\ 

That is, F XiXj (AP ) < 0 precisely when b G [(1 + k)|P|, \P\\. 

On the other hand, if A G (1/(1 + k), +oo), then A > 1 since -1 < k < 0. If AP is 
on the surface, which means A|P|+k (A - 1)|P|= b, we have 

b + k |P| 


A 


(1 + k) |P| 


So A > 


1 


1 + K 


implies 


b > (1 - k)|P|. 

But (1 - k)|P|> |P| and by the physical restriction on b, we must have k|P|< b < \P\. 
So we see that 

\P\< (1 - k)|P|< b < |P| 

which is impossible. 

b + k |P| 


Let now A G (1,1/(1 + k)). Then since A 


(1 + k) |P| 


we see that 


, b + k |P| 1 

1 < --“7T7 < 


(1 + k) |P| 1 + K 


so that 


\P\< b < (1 - k)|P| 

But by the constraint on b we must have b < |P|, which contradicts the above 
inequality 

Hence given the critical values of b as discussed above, the main result in this 
section is as follows: 
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Theorem 3.4. Suppose -1 < k < 0. The surface |X|+k |X - P |= b is 

(1) convex, when b e [k\P\, (1 + k)|P|]; 

(2) neither convex nor concave, when b £ ((1 + k)\P\, |P|). 

Ifx<-1 then the surface |X|+k |X - P|= b is 

(1) convex when b e ((1 + k)\P\, |P|); 

(2) neither concave nor convex when b e (k\P\, (1 + k)\P\], 


Proof We begin by proving Theorem 3.4 in dimension two when — 1 < k < 0. We 
will determine the curvature using Proposition 3.1]. That is, we let £ = (-F y , F x ) 
into the quadratic form Q(X, £) and determine its sign for each X on the curve 
when b varies. We have that 


Q(X,(-F y ,F x )) = | X _ p|3|x|3 X 
[k 3 |X| 3 |X-P| 2 +|X| 2 |X-P| 3 

+k 2 |X - P|(X • P - 3|X| 2 ) (X • P - |X| 2 ) + k|X|(|X| 2 -X • p) (2|P| 2 -5X • P + 3|X| 2 )] 

: = |x _p |3|x|3 x [k 3 |X| 3 |X - P| 2 +|X| 2 |X - P| 3 +B + C]. 

Since 2(X • P) = |X| 2 +|P| 2 -|X - P| 2 , we can write B, C as functions of |X| and |X - P| 
as follows: 


B = k 2 \X - P|(3|X| 2 -X • P) (|X| 2 -X • p) 

= k 2 |X - P|(|X| 2 -X • pf + 2k 2 |X - P| |X | 2 (|X| 2 -X • P) 

= i k 2 |X - P|(|X| 2 +|X - P| 2 -|P | 2 ) 2 + k 2 |X - P\ |X | 2 (|X| 2 +|X - P| 2 -|P| 2 ); 


C = k|X|(|X| 2 -X • P) (2|X - P| 2 +|X| 2 -X ■ P) 

= 2 k|X| |X - P| 2 (|X| 2 -X • P) + k|X|(|X| 2 -X • pf 
= k|X| |X - P| 2 (|X| 2 +|X - P| 2 -|P| 2 ) + i k|X|(|X| 2 +|X - P| 2 -|P| 2 ) 2 . 
The numerator of Q(X, (—F y , F x )) then equals 


|X| 2 |X - P| 2 (k 3 |X|+|X - P\) + (|X - P| 2 +|X| 2 -|P | 2 ) 2 (|X|+k|X - P|) 

+ k|X||X - P\(\X - P| 2 +|X| 2 -|P| 2 ) (k|X|+|X - P|) := H(|X|, |X - P|, k, \P\), 
and we want to optimize this quantity when X is on the curve F(X) = |X|+k|X-P|= 
b. If X is on the curve there holds |X - P|= --—- > 0. Hence |X|> b, and also 


K 
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||P|-|X|| < |X- P|= ^—^3. Hence 


K 


6 - 1X1 < |P|-|X|< 6 - 1X1 


K 


K 


Multiplying by k < 0 yields 


so 


and 


-& + |X|> k\P\-k\X\> b - |X|, 
(1 - k)|X|> b - k\P\, 


(1 + k)|X|> b + k\P\. 

Then the last two inequalities together with |X|> b imply that X must satisfy 

b - k\P\ b + k\P\ 


|X|> max • 


1 — K 1 + K 


-M- 


In addition, from |X - P\= --—- we get that |X|+|P|> |X - P|= --—-, which 

K 

implies 


K 


|X|< 


b - k|P| 


since -1 < K < 0. Substituting |X - P\= 


1 + K 

b-\X 


K 


in the expression for H above we 


obtain the following expression in terms of K, b, |X| and \P\: 
b{b 2 - k 2 |P| 2 ) 2 + 6 b{x 2 - 1 )(b - x\P\)(b + k\P\)\X\ 2 

4 k 3 

4(k 2 - 1) ( b 2 (k 2 - 2) + k 2 |P| 2 ) |X| 3 -3&(k 2 - 1) 2 |X | 4 v| 7 x 
+- 4^3 - : = G(\X\,b, k), 

with |X|> b. Thus the problem of finding the b's for which Q(X,(-F y ,F x )) > 0 
for all X, reduces to minimizing a polynomial of one variable |X| over the range 


max 


b - k\P\ b + k\P\ 


1 — K 1 + K 

b - k\P\ 


,b\< |X|< 


b - k\P\ 
1 + k 


. Notice that the value of the maximum 


is given by M = 


1 — K 


The derivative of G with respect to the first variable 


b — k\P\ 

This means that the ball centered at zero with radius —-contains the oval. In fact, this is the 


1 + K 

smallest ball centered at zero that contains the oval because the point X = AP with A = 


b - k|P| 
(1 + k)|P| 


-b - \P\ 

satisfies the equation F(X) = b (notice that A < 0 and A - 1 = —-—— < 0 since —b < -k|P|< |P| 

(1 + k)|P| 

since —1 < k < 0). 

"Here we use that b < |P|, because the point P must be outside the oval to have a problem 
physically meaningful. This means that the radius of the largest ball centered at zero and contained 
b - k|P| 


in the oval has radius 


1 - K 
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t = |X| equals 


d t G(f, b, k) 


3(-l + K 2 )(b - t)t(b 2 - K 2 \P\ 2 +b(-l + k 2 )0 

¥ 


_^2 + K 2 | p |2 

and therefore if b # 0 then c) t G{t, b, k) = 0 when t = 0, b, or p := . If 

b(-1 + k 2 ) 

b = 0 then d t G(t,0,K ) = 0 when t = 0. Furthermore, the coefficient of t 4 in the 

3b 

original function G(t,b,k ) equals L b := -- —-(1 - k 2 ) 2 . 

4 K 

To optimize G(t, b, k ) we proceed by analyzing the following cases. 


Case k\P\< b < 0. 

Since k < 0, then the coefficient L b < 0. We also have t b > 0. We write 
d t G(t, b, k) =-3 b (K ' 2 ~ 1)2 t(t-b)(t- t h ). 

For t b < t < co we have d t G(t, b, k ) < 0, and for 0 < t < t b we have d t G(t, b, k) > 0. 
We optimize the function G on the interval 


(b - k\P\ b-x\P\\ 
\ 1-K ' 1 + K ) 


C (0, tb). 


We have that d t G(t, b, k) > 0 for tel and so G is strictly increasing on I. Hence 
we obtain 


max G(t, b,K) = G 

tel 


b - k\P\ 

1 + K 



(b + \P\) 2 (b - K\P\) 2 (b + \P\-k\P\) 

(T TW 


and 


min G(t, b,x) = G 

tel 


b - k\P\ 

1-K 



(b - \P\) 2 (b - x\P\) 2 (-b + (1 + k)\P\) 


Note that this minimum is positive since -1 < k < 0 and the maximum is 
positive as well since b > k|P|. Hence the maximum and minimum have the same 
sign and as such the curve is convex for this range of b. 


Case 0 < b < (1 + k)|P|. 

By calculation we have that L b > 0. Let us analyze the sign of t b . Since b > 0 
and -1 < k < 0, we have tb > 0 if and only if b > — k|P|, and t b < 0 if and only if 
b < —k|P|. Let us compare — k|P| with (1 + k)|P|. Wehave-x'|P|< (l + K)|P|ifandonly 
if k < -1/2, and -k|P|> (1 + k)|P| if and only if k > -1/2. Therefore, if k > -1/2 
then t b < 0. On the other hand, if k < -1/2, then b e (0, — k|P|) U [—k|P|, (1 + k)|P|). 
If b e (0, —k|P|), then t b < 0. If & 6 [—k|P|, (1 + k)|P|), then t b > 0. Notice also that 
t b < b since b < (1 + k)|P|. If k > -1/2 then t b < 0 and d t G < 0 on (0, b) while d t G > 0 
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on (b, +oo). Since 


b - x\P\ 


x 


> b, we have that the interval I C ( b , +oo), and 


min G(t, b,x) = G 

tel 


'b - k\P\ 

, 1 - K 


,b,K 


(i b - \P\) 2 (b - x\P\) 2 (-b + (1 + k)\P\) 


(-1 + kf 

If K < -1/2 then 0 < tj, < b and so d t G > 0 on (0, tj,), d t G < 0 on (fj,, b) and d t G > 0 
on (b, +oo). Once again 1 C (b, +oo), and we obtain 


min G(f, b,x) = G 


'b - k\P\ 


tel 


K 


,b,K = 


(b - \P\) 2 (b - x\P\) 2 (-b + (1 + k)\P\) 

FTTf ' 

Note that this minimum is positive since b < 0. Furthermore we have 

\ r^( b ~ K l P l 7 \ 

max G (t, b.x) = G —-, b, k 

tel \ 1 + K ) 

which is positive because b > k|P|. Hence the maximum and minimum have the 
same sign and for this range of b' s the curve is also convex. 


Case (1 + k)\P\< b < \P\. 

In this case we have that 4 < b < 
'b-x\P\ b - k\P\ 


b - k\P\ b - k\P\ 


so that the interval I 


1 - k 1 + K 

C ( b,oo ). Furthermore we have that 4 > 0 provided that 


\ 1 - k 1 + K 

—1 < K < —1/2 and -k|P|< b < |P| or if -1/2 < k < 0 and (1 + k)|P|< b < |P|. 
Moreover we have that 4 < 0 if —1 < k < —1/2 and (1 + k)|P|< b < —k|P|. 

Also on the interval 1 we have that d f G(t, b, k) > 0 so that 


min G(f, b, k) = G 

tel 


'b - k|P| 


k 


,bx = 


(b - \P\) 2 (b - K\P\) 2 (-b + (1 + k)\P\) 
(-1 + k ) 2 


which is negative since b > (1 + k)|P|. Furthermore we have that 


max G(f, b,x) = G 


tel 


'b-x\P | 
i 1 + x 


,b,x = 


(b + |P|) 2 (b - K|P|) 2 (b + (1 - x)\P\) 
(1 + k ) 2 


which is positive since -1 < k < 0. Thus the minimum and maximum of G on the 
interval have opposite signs which implies that the curve is neither convex nor 
concave for this range of b. 

In dimension three, the above result still holds true because the oval is radially 
symmetric with respect to the axis OP. 

Finally suppose that x < -1. 

We have X £ {X : |X|+k|X-P|= b} := Oif andonlyifX-P e jz : |Z|+^|Z + P|= ^ j 
O'. The ovals O and O' have the same curvature, and by the -1 < x < 0 case. O' is 
convex provided — |P|< — < fl + — ] |P| and neither convex nor concave provided 

K K \ 1C / 
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(l + — W|< — < \P\. Multiplying both inequalities by k we obtain the desired 

\ K1 K 

range of b. 

□ 


4. Fresnel Formulas for NIMs 


To obtain the Fresnel formulas for NIMs materials we briefly review the cal¬ 
culations leading to the Fresnel formulas for standard materials, see Q1 Section 
1.5.2], 

The electric field is denoted by E and the magnetic field by H. These are three- 
dimensional vector fields: E = E(r, t ) and H = H(r, f), where r represents a point in 
3 dimensional space r = ( x , y,z). The way in which E and H interact is described 
by Maxwell's equations QJ and |[12l Sec. 4.8]: 


(4.28) 

(4.29) 

(4.30) 

(4.31) 


VxE = --§' 

tt 271 edP 

V X H = —crE + - — 

C C dt 

V • (cE) = 47ip 

V • (pH) = 0, 


c being the speed of light in vaccum. 

We consider plane waves solutions to the Maxwell equations (4.28), ( 4.29| ), 
(4.30|), (4.31) with p = 0 and o = 0, and having components of the form 


r • s 


a cos \co[t -+ 5 = a cos (cot - k • r + 6) 


co 


with k = —s and a, b real numbers, and s is a unit vector. If E(r, t) = E(k ■ r - cot ) 

and H(r, t ) = H(k • r - cot) solve the Maxwell equations ( 4.28| ), ( |4.29| ), ( 4.30 ), one 
can show that 


(4.32) 


H = — (k x E), and , E = (k x H). 
pco eco 


We assume that the incident vector s = s 1 with 

A. /V 

s = sin @,i + cos @,k. 


That is, s' lives on the xz-plane and so the direction of propagation is perpendicular 
to the y-axis. Also the boundary between the two media is the xy-plane, v denotes 
the normal vector to the boundary at the point P, that is, v is on the z-axis, and 
is the angle between the normal vector v and the incident direction s' (as usual 
i,j,k denote the unit coordinate vectors). Let I\\ and I ± denote the parallel and 
perpendicular components, respectively, of the incident field. The electric field 
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corresponding to this incident field is 
(4.33) E'(r, t ) = (—J|| cos 6j,I ±/ I\\ sin 0 Z ) cos 
with 


r • s 

Vi 


= E^ cos [t 


r • s 


Vi 




From ( 4.321 ), the magnetic field is then 

H'(r,0= a &xE ! ) 

V Fi 


(-Jj_ cos 9i, -I\\, I ± sin 0,) cos Ico It - 



r • s 

V\ 


Hq cos [co \t - 


r • s 

Vi 


Let us now introduce s f , the direction of propagation of the transmitted wave, 
and let 0 f be the angle between the normal v and s f . Similarly, s' is the direction 
of propagation of the reflected wave and 0,. is the angle between the normal v and 
s''. We have that s' = sin 0, i + cos 0, k = sin 0,i - cos 0,k. Then the corresponding 
electric and magnetic fields corresponding to transmission are 


E f (r, t) = (-Ty cos 0 t , T ±r Ty sin 9 t ) cos \co it - 


r ■ s 


v 2 


En cos [colt - 


r • s 


v 2 


H f (r, t) 
with 



— (-Tj_ cos 0 f , —Tn, T ± sin 0 f ) cos \co f - 


r ■ s 


v 2 


H,( cos [co[t - 


r ■ s 


v 2 


v 2 


yfeiiii 


There are similar formulas for the fields E' (r, f), H' (r, t) corresponding to reflection. 
Since the tangential components of the electric and magnetic fields are continuous 
across the boundary, see e.g. HI Section 1.1.3], we obtain the equations 


I± + R± - Tj_, 


cos 9i(I\\ - R\\) = cos 0fT||; 


and 


In R\\ 

+ 


Tn 





COS 0; 


The wave impedance of the medium is defined by 


R i 




= cos 9 t 


T, 



z = 
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r / 

—, j = 1,2, and solving the last two sets of linear equations 


(4.35) 


Rn = 


2 Z\ cos Gi j 

Z2 COS Gi + Zi cos G t 11 

z 2 cos Gj - Zi cos G t i 

z 2 cos Gi + Zi cos G t 11 


Ti 

R 


2 Zi cos 0, 

Zi cos Gi + z 2 cos G t ± 
Zi cos Gj - z 2 cos 0 f 
Zi cos 0, + z 2 cos G t 1 


These are the Fresnel equations expressing the amplitudes of the reflected and 
transmitted waves in terms of the amplitude of the incident wave. 

We now apply this calculation to deal with NIMS. We will replace s 1 by x and 
s* by m, and we also set 

V £ 2/-l2 

K =- 

Recall v is the normal to the interface. We have cos 6, = x ■ v and cos 6 t — m ■ v. 
Suppose medium I has €\ > 0 and pi > 0; and medium II has e 2 < 0 and p 2 < 0. In 
other words, medium I is "right-handed" and medium II is "left-handed". This 
means in medium II the refracted ray m stays on the same side of the incident ray 
x. 


From the Snell law formulated in (|2.7| (notice that k < 0), we have x - Kin = Av, 
so the Fresnel equations (4.34|> and (4.35) take the form 


(4.36) 


T II 

T ± 

R\\ 

Rx 


2 Zi x • v 

(z 2 x + zi m) ■ 

• V 

2 zi x ■ v 


(zix + z 2 m) ■ 

■ V 

(z 2 x - Zi m) ■ 

■ V 

(z 2 x + zi m) ■ 

■ V 

(zix - z 2 m) ■ 

■ V 

(zix + z 2 m) ■ 

■ V 


2 z\ x ■ (x - k m) i 
11 (z 2 x + Z\in) • (x — k m) 11 
2 z\ x ■ {x - k m) j 
± (zix + z 2 m) • (x - k m) 

I (z 2 x -Z\m) ■ (x- k m) 

11 (z 2 x + Z\in) ■ (x — k ni) 11 
{z\X - z 2 m) -(x-Ktn) 

± (zix + z 2 m) ■ (x - k m) 


Notice that the denominators of the perpendicular components are the same and 
likewise for the parallel components. 

As an example suppose that e\, pi are both positive and e 2 = -e 1 and p 2 = -px. 
This is the so called mirror like material. Then Zi = z 2 and k = -1, and so 
R|l = R ± = 0 for all incident rays. This means that all the energy is transmitted 
and nothing internally reflected. Notice also that if k / -1 then there is a reflected 
wave, that is, R ± and R\\ may be different from zero. 


Remark 4.1. The wave impedance, like the index of refraction, is a unique char¬ 
acteristic of the medium in consideration. However, unlike the refractive index n, 
the wave impedance z remains positive for negative values of e l , p ; . Also, compare 
equations ( 4.35| ) with the table in Ifl5l pg. 765] giving R ± . There the Fresnel for¬ 
mula for r ± in the exact Fresnel formula column contains a misprint. We believe 
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our calculations above to be correct and consistent with what is expected in the 
case of negative refraction. 

4.1. The Brewster angle. This is the case when R\\ = 0, this means when 

Z 2 cos 9i - Z\ cos 6 1 = 0, 


which together with the Snell law (sin 6 , = k sin 6 t ) yields 


—) cos 2 9j + ~r sin 2 9j = 1 = cos 2 0, + sin 2 9, 

Z\) K z 


So 


tan 2 9i = 


U -i 

Zl 


1 - 


_ _ /A / £lR2 ~ g2/A 

l_ H 1 \e 2 p 2 -GijUy 

K 2 


if 1 /k 2 # 1, and the Brewster angle is 0, with 


f- 12 1 

( UP2 - e 2j Ui \ 

1 1-9 ' 

ie 2 p 2 -eipi/ 


Let us compare this value of 0, with the case when k > 0. As above the Brewster 
angle occurs when R\\ - 0. In this situation, the reflected and transmitted rays are 
orthogonal to each other, and via Snell's law, it follows that 


tan 0/ = n 2 /fti 


which is consistent with what we would expect. 


4.2. Fresnel coefficients for NIMs. To calculate these coefficients we follow the 
calculations from flv Section 1.5.3]. The Poynting vector is given by S = —ExH, 

where c is the speed of light in free space. From ( 4.32| > we get that 

fe". 


—Ex(—kxE) 

471 jJ.0) 


c 

471 


-E x s x E, 


co 


where k = — s, s a unit vector, and v = —— 

V yjefl 


. Using the form of the incident wave 


( 4.33| ), the amount of energy /' of the incident wave E' flowing through a unit area 
of the boundary per second, is then 


f = IS'lcos 0, 



Similarly, the amount of energy in the reflected and transmitted waves (also given 
in the previous section) leaving a unit area of the boundary per second is given 
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by 


/ r = is'icos = AL /—|e, 

471 \ fl , 
^2, 


,q| 2 cos Q l 


/ f = |S f |cos@ f =— /-|E'| 2 cos0 f . 
471 Y H2 

The reflection and transmission coefficients are defined by 


* = £ = (!5!' 

V VIE'U 


P 

, and T = — 

V 


je 2 /ti cos 0 t /|Eq| 
Ci/.i 2 cos 0,- \|EU 


By conservation of energy or by direct verification, we have ( R + 7~ = 1. 

For the case when no polarization is assumed, we have from Fresnel's equations 
( 4.36| > that 


|E 0 i - 

and so 
<R = 


n2 ~ R 2 + R 2 


'[5 

JE' 


(z 2 x — Z\ vi) ■ (x — k m ) 
(z 2 x + Z\ m) ■ (x - k m) 

2 R; + R 2 


-.2 


In + 


(zi x - z 2 m) ■ (x - k m) 
(zi x + z 2 m ) ■ (x - k m) 


A. 


f 2 + f 2 


(z 2 x - Z\ m) ■ (x - k m) 
(z 2 x + zi m) ■ (x - k m) 


f 2 


f 2 + f 2 


+ 


(zi x - z 2 m) ■ (x - k m) 
(zi x + z 2 m) • (x - k m) 
2 


f 2 


f 2 +f 2 



which is a function only of x • m. We assume in these formulas that the product 
c,/i, > 0, / = 1,2. In principle the coefficients Zy and Z ± might depend on the 
direction x, in other words, for each direction x we would have a wave that changes 
its amplitude with the direction of propagation. The energy of the incident wave 
would be f(x) = |E^| 2 = I||(x) 2 + Z ± (x) 2 . Notice that if the incidence is normal, that 
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(Z 2 — Z\ y 

is, x = m, then R = - which shows that even for radiation normal to the 

\z 1 +z 2 ) 

interface we may lose energy by reflection. 

Remark 4.2. The Fresnel coefficients R and T being written in terms of x and m 
as above are useful for studying some refraction problems in geometric optics. 
In fact, when k > 0, it was shown in ||6J that it is possible to construct a surface 
interface between media I and II in such a way that both the incident radiation 
and transmitted energy are prescribed, and that loss of energy due to internal 
reflection across the interface is taken into account. 

5. Conclusion 

We have given a vector formulation of the Snell law for waves passing between 
two homogeneous and isotropic materials when one of them has negative refrac¬ 
tive index, that is, is left-handed. This formulation was used to find surfaces 
having the uniform refraction property both in the far field and near field cases. 
In the near field case, in contrast with the case when both materials are standard, 
these surfaces can be neither convex nor concave and can wrap around the target. 
A quantitative analysis in terms of the parameters defining the surfaces has been 
carried out. We have used the vector formulation of Snell's law to find expressions 
for the Fresnel formulas and coefficients for NIMs. We expect these formulas to 
be useful in the design of surfaces separating materials, one of them a NIM, that 
refract radiation with prescribed amounts of energy given in advance. This will 
be done in future work and requires more sophisticated mathematical tools. 
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